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. Abstract 

m 



A reduction of Benney's equations is constructed corresponding 



ly-) ■ to Schwartz-Christoffel maps associated with a family of genus six 

, cyclic tetragonal curves. The mapping function, a second kind Abelian 

' integral on the associated Riemann surface, is constructed explicitly 

as a rational expression in derivatives of the Kleinian er-function of 
the curve. 

> ■ 

X 

^ : 1 Introduction 

In [1, 2] it was shown that Benney's equations 

A n t = A n x +l + nA n ~ l A x} n > 0, 

admit reductions in which only finitely many N of the moments A n are 
independent, and that a large class of such reductions may be parametrised 
by conformal maps from the upper half p-plane to a slit domain - the upper 
half A-plane, cut along N nonintersecting Jordan arcs, which have one fixed 
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end point on 9f(A) = 0, and whose other 'free' end is a Riemann invariant of 
the reduced equations. 

A natural subclass of these occurs where these Jordan arcs are straight 
lines, leading to a polygonal domain and hence an A-parameter Schwartz- 
Christoffel map; an important and tractable subfamily of these is the case in 
which the angles are all rational multiples of tt, and in this case the mapping is 
given by an integral of a second kind Abelian differential [3] on an algebraic 
curve. Such examples have been worked out explicitly, in [3], [4], [5] and 
[6]. These have looked at elliptic and hyperelliptic curves as well as a cyclic 
trigonal example. It is thus worthwhile to generalise this to other algebraic 
curves. 

In all these examples, reductions have been constructed explicitly - both 
the integrand and its integral were evaluated using quotients of derivatives 
of the a-function associated with the respective curves. These curves are 
all specific examples from the wider class of cyclic (n, s) curves which have 
equations of the form: 

y n = x S + fls-iX 3 ' 1 H VUlX + flQ. (1) 

We suppose that (n, s) are coprime with n < s, in which case the curves have 
genus g — |(n — l)(s — 1), and a unique branch point oo at infinity. In this 
paper we consider a reduction associated with a cyclic tetragonal curve, that 
is, from the class (1) with n = 4. For simplicity we look at the case with 
s = 5 here. 



2 Benney's Equations 



In 1973, Benney considered an approximation for the two-dimensional 
equations of motion of an incompressible perfect fluid under a gravitational 
force [7]. He showed that if moments are defined by: 



A n (x,t) = [ u n dy, 
Jo 



where u(x,y,t) is the horizontal fluid velocity and h(x,t) the height of the 
free surface, the moments A n (x,t) satisfy an infinite set of hydrodynamic 
type equations 

dA n dA n+1 dA 

^r + — +nAn - l ^ x - = (» = 1 ' 2 .-), (2) 
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now called the Benney moment equations. 

Identical moment equations can alternatively be derived from a Vlasov 
equation [8], [9]: 

^ + p<V_^W = o (3) 

dt 2 dx dx dp 

Here / = f(x,p,t) is a distribution function and the moments are defined 
instead by 

poo 

A n = / p n fdp. 



We assume throughout that / is such that all these moments exist. The 
equation of motion (3) has the Lie-Poisson structure : 

d -±Ai 5 4\ =<>. ( 4) 



at V'sfj 



p,x 



where {•,•} is the canonical Poisson bracket. Kupershmidt and Manin 
showed directly that the moment equations are Hamiltonian [10], [11]. If we 
set H =\H 2 = \(A 2 + A 2 ), A = (A , A u . . .), then 

d A- B d JL (5) 
dt~ B dA {b) 

where the matrix operator B is given by 

d 8 

£>n,m — n A n + m ^i— h m—— ■ A n + m -i. 

OX OX 

This is consistent with (4) in the sense that if H is some function only of the 
moments, the moment equations resulting from (4) and (5) are identical. 

Benney showed in [7] that system (2) has infinitely many conserved den- 
sities, polynomial in the A n . One of the most direct ways to calculate these 
is to use generating functions [10]. Let X(x,p,t), a formal series in p, be the 
generating function of the moments 

A(*,p,t)=P + E^r W 

n=0 y 

and let p(x, X,t) be the inverse series 



OO 

p(x,X,t) A 



m=0 
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We note here that if A n = f^° 00 p n fdp is substituted into (6), then this can 
be understood as the asymptotic series, as p — > oo, of an integral 

Here p' runs along the real axis, and we take $s(p) > 0. It follows that 
X(p) is analytic in its domain of definition. If f{p) is Holder continuous, the 
boundary value, on the real p-axis, of \{p) will itself be Holder continuous. 
Comparing the first derivatives of X(x,p,t), we obtain the PDE 

d X dX dX ( dp dp dA \ 
dt ^ dx dp \dt ^ dx dx ) ' 

If we now hold p constant, this gives 

5A +p 5A_^4o5A =Q (g) 
dt dx dx dp 

which is a Vlasov equation of the same form as (3). Thus (3) and (9) have 
the same characteristics. Any function of A and / must satisfy the same 
equation. 

Alternatively, if we hold A constant in (8), then we obtain the conservation 
equation 

dp d ( 1 2 



Substituting the formal series of p(x,X,t) into (10), we see that each H n is 
polynomial in the A n and is a conserved density. Any of the H n could be 
used as the Hamiltonian in (4), and the resulting flows all commute. From 
this we define the Benney hierarchy to be the family of evolution equations 

dt n { n 6f J 
Again, A satisfies an equation analogous to this, 
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so that both / and A are advected along the same characteristics. These 
characteristics are flows of a Hamiltonian vector field, with Hamiltonians 
given by the relation: 

l5H n \ = /An 

n Sf J \n) + 
where (•)+ denotes the polynomial part of the Laurent expansion. 

2.1 Reductions of the moment equations 

Suppose that for some family of points, p = Pi(x, t), X(pi) = A«(x, t), we have 

dX 
dp 



p=Pi 



Then (8) reduces to: 



dt % dx 



where %r = %r - and = . . We say that A; is a Riemann invariant 

at at I p=pi ox ox I p=pi J 

with characteristic speed pi. We will see that there are families of functions 
X(p) which are invariant under the Benney dynamics, and are parametrised 
by N Riemann invariants Aj. 

A hydrodynamic type system with N > 3, independent variables can not 
in general be expressed in terms of Riemann invariants. If such a system 
does have N Riemann invariants, it is called diagonalisable. Tsarev showed 
in [12] that if a diagonal hydrodynamic-type system 

^ + ^(A)^ = (, = 1,2,..., AT). (11) 

is semi-Hamiltonian, that is if 

diV k \ n / djV k 



for k distinct, where 



d k = . 
dX k 
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then it can be solved by the hodograph transformation. Any Hamiltonian 
system of hydrodynamic type is semi-Hamiltonian. Given a second equation 
of type (11) 

^+^(A)^ = (i — 1,2, ... , N), (12) 

and requiring it to be consistent with (11), we find that the Wi(X) must satisfy 
the over-determined linear system 

d k Wi d k Vi 

l f: k. (13) 



Wi -w k Vi- v k 



These equations are consistent provided (11) is semi-Hamiltonian. If the 
condition (13) holds, we say that (11) and (12) commute. In this case a set 
of equations for the unknowns \(x,t) is given by : 

Wi(\)=Vi(\)t + x, (i = l,2,...,N) 

where t and x are the independent variables. Thus any reduction of this type 
can be solved in principle. 

This generalized hodograph construction cannot easily be applied directly 
to the Benney equations however, as these have infinitely many dependent 
variables. Instead we will now consider families of distribution functions /, 
which are parameterised by finitely many N Riemann invariants \(x, t). We 
are interested in the case [1], [2] where the function X(p,x,t) is such that 
only N of the moments are independent. Then there are N characteristic 
speeds, assumed real and distinct, and N corresponding Riemann invariants 
(pi, Aj), so Benney's equations reduce to a diagonal system of hydrodynamic 
type with finitely many dependent variables Aj, 

w +mX) i£ = ( < = 1 » 2 »-» jv )- ( 14 ) 

Such a system is called a reduction of Benney's equations. 

The construction of a more general family of solutions for equations of 
this type was outlined in [1] and [2]. An elementary example is the case 
where the map A + takes the upper half p-plane to the upper half A-plane 
with a vertical slit as follows. This is a Schwarz-Christoffel map: 

p f — Pi 

X + (x,p,t)=p+ ==dp. 

Joo V(p '-Pl){p '-P2) 
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If the residue at infinity is set to be zero, then this imposes the condition 
Pi — \(Pi +P2) an d we get the solution 



\+(x,p,t) = Pi + y/p 2 - (pi + P2)P + PlP2 

= Pi + \ftp~ ~ Pi) 2 + 2A 

(from the expansion as p — > 00). This gives a steadily translating solution of 
Benney's equations (4) 



dt 



p,x 



The two parameters p\ and p 2 are not independent, as for consistency their 
sum must be a constant. Hence only the end point of the slit in the A-plane 
is variable. This is the Riemann invariant. 

This construction can be generalized, [2], mapping the upper half p-plane 
to the upper half A-plane with N curvilinear slits. However, in this paper we 
are specifically interested in straight slits - these mappings are all of Schwartz- 
Christoffel type. 



2.2 Schwartz- Christoffel reductions 

The case of a polygonal TV-slit domain is of particular interest. The real 
p-axis has M vertices Uj marked on it; the preimage in the p-plane of each 
slit runs from a vertex pj, to a point Vi, the preimage of the end of the slit, 
and then to another vertex pj+i- The angle 7r in the p-plane at pj is mapped 
to an angle ajir at the image point. The internal angle at the end of each 
slit is 2n. 

The mapping function is then given, up to a constant of integration, by 



A 



Ui=i(p-^i 

UZ(p-PjY 



dp. 



If the integrand is to converge to 1 as p — > 00, we require 



to avoid a logarithmic singularity, we further impose Y^jLi a jPj = 
We then define A more precisely as 



2JV 
3 



1 a 3 



N, while 



N 



A 



/v 
-c 



Uf=iiP-Vi) 

nJliCp-ft-) 1 - ' 



1 



dp. 
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Other constraints are imposed by requiring the vertices pj to map to points 
A°, the fixed base points of the slits; there remain N independent parameters, 
which can be taken to be the movable end points of the slits Xi(x, t). These 
satisfy the equations of motion 

dt 1 dx 

To understand and to solve these equations, it is necessary to understand the 
dependence of the Vi on the Riemann invariants (Ai, . . . , A at) - we thus need 
to evaluate the Schwartz- Christoffel integral explicitly. 

The most tractable cases are where all the ctj are rational, so that the 
integrand becomes a meromorphic second kind differential on some algebraic 
curve. In this case the only singularity is as p — > oo, where the integrand 
has a double pole with no residue, and the integral thus has a simple pole. 
For specific families of curves, as in [3], [4], [5] and [6], this integral has been 
worked out explicitly. In each case, these mappings were found as rational 
functions of derivatives of the Kleinian a-function of the associated curve. 



3 A tetragonal reduction 

We will consider reductions that allow us to work on a tetragonal surface; 
these have not been considered before. Such reductions will require two or 
more sets of straight slits, making angles of 7r/4, tc/2, 3tc/A to the horizontal. 
Define V to be the upper half p-plane with 14 points marked on the real axis, 
as in Figure 1. 



4 ■ »- 




Figure 1: The domain V within the p-plane. 



These points satisfy 

pi < Vi < p 2 < V 2 < p3 < h < p A < p 5 < V 4 < p6 < V 5 < p 7 < V 6 < p 8 - 

Then define the domain C as the upper-half A-plane with two triplets of slits, 
as described above. We let the first trio of slits radiate from the fixed point 



8 



Figure 2: The domain £' within the A-plane. 



Pi, with the end points of these three slits labelled Vi,v 2 and v s respectively. 
Similarly, let the second trio of slits radiate from p§ and have end points 
^4,^5,^6- Finally impose the conditions that 

A(p) =p, Pl =p 2 =p 3 = P4 , 

\(y) =V, P5=P6=P7=P8- 

We then see that C! is the slit domain as shown in Figure 2, and the mapping 
A : V — > C! can be given in Schwartz Christoffel form by 

\(p)=p+ f [yip') - l]dp' (15) 

J oo 



where 



where 



<p(p) = n? =1 (P-^) = yiUp-v.) 



V A = tl(p-iH). (17) 

i=i 

Note that we require the following zero residue property: 

]wKp(p)~l + o(\). (18) 

p^oo \P J 



This mapping would lead us to consider the Riemann surface given by 
points (p,y) that satisfy (17). However, we wish to consider the simplest 
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Figure 3: The domain C within the A-plane. 



possible tetragonal surface (one with only six branch points) and so we col- 
lapse two of the slits, (the final two by choice). This simplifies our A-plane 
to C, given in Figure 3. 

Further, as in the trigonal case, the analysis of this surface is eased if 
we put it into canonical form, by mapping one of the branch points (p 8 by 
choice), to infinity. So we use the following invertible rational map to perform 
these simplifications on our curve and integrand. 



Pa=P8, P7=Ps, v 5 =p s , v 6 =ps, 
p = p 8 -(!/*)» Pi=P&-w., i = l,..., 5 



y 



sk 



5 5 ^ 

where k A = - JJ(j3 8 - p t ) = - J[ — 



i=i i=i 
If we perform the mapping (19) on the curve (17) we obtain 



(19) 



s A k A 

IF 



IJ ( P8 - 7 - PS + 

.1=1 



t 



Ti 



1 



Ps-j-Ps 



t[(t-Tj 

.i=i 



•(-1) 



5 1 

.i=i 



This simplifies to give the following canonical form of (17). 



s 4 = n (t - =t 5 + /i 4 t 4 + fi 3 t 3 + /i 2 t 2 + ntf + no, 



(20) 



i=i 
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for constants /i , ...,114. Let C denote the Riemann surface defined by (20). 

We now consider A(p) as mapping V — > C by performing (19) on the 
integrand (16). 




= K [A 4 t 4 + A 3 t 3 + A 2 t 2 + ^t+l]~E (21) 

where = — A/k 3 and Ai, . . . , A4 are constants. We will evaluate this inte- 
grand using algebraic functions defined upon C. 

4 Properties of the tetragonal surface C 

The Riemann surface C is defined by (20), the cyclic tetragonal curve of 
genus six. This is also referred to as the cyclic (4,5)-curve. The surface is 
constructed from four sheets of the complex plane, with branch points of 
order 4 at 7\, . . . , T 5 , T 6 — 00, a local coordinate at t = 00 given by £ = t' 1 ^ 
and branch cuts along the intervals 

[T U T 2 ], [T 2 ,T 3 ], [T 3 ,T 4 ], [T 5 ,oo]. 

This surface was recently considered in [13] , where the aim was to generalise 
the theory of the Weierstrass p-function to Abelian functions associated with 
(20). We will give the essential properties of the surface here, but refer the 
reader to [13] for some of the details and proofs. 

We start by noting that there are a set of Sato weights associated to 
the surface, which render all equations within the theory homogeneous. The 
weights of the variables t, s and the curve constants can be determined up 
to a constant factor, from the curve equation. We may set the weights as 
below. 







s 


fJL 4 


A*3 


£*2 




IM) 


Weight 


-4 


-5 


-4 


-8 


-12 


-16 


-20 
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The weights of other variables and function in the theory are derived uniquely 
from these, with all other constants are assigned zero weight. 
We define a basis of holomorphic differentials upon C by 



du = (dui, . . . , cIuq), duAt, s) = s ) dt, 



9i(t,s) = 1, 

04(*,S) = t 2 , 



4s a 
g 2 (t, s) = t, 
g 5 (t,s) = ts, 



9s(t, s) 
9a(t: s) 



where 

We can use the local parameter £ to express these as series, 



s, 

„2 



du 1 = -e°+o(e i )dti 
du 2 = ~e+o(en 
du 3 = -e+o(e)dt 



du A = -e + \^e+o{e)di 
du 5 = -e + \^e+o{e)di 
du 6 = -i+\t M e+o(e)dc 



(22) 



(23) 



We know from the general theory that any point u e C 6 can be expressed as 



U = («i, U 2 , U 3 , W 4 , Ub, U 6 ) = / dU ' 

i=l 



where the Pi are six variable points upon C . Integrating (23) gives 



«i = "ire 1 + 0(£ 15 ) u 3 

u 2 = - i 1 e+o(e i ) u 4 



-le+o(e) u 6 



L 2 e+o(e 



(24) 



from which we can conclude that the weights of u are: 





Ui 


u 2 


u 3 


«4 


W 5 


Ml 


Weight 


+ 11 


+7 


+6 


+3 


+2 


+ 1 



Next we choose a basis of cycles (closed paths) upon the surface defined by 
C. We denote them 

&i,Pj, 1 < i, J < 6, 
and ensure they have intersection numbers 



di ■ OLj = 0, Pi • Pj = 0, 



OL-l 



pj = Sij = I 



1 if i = j 

o if i ^ j ■ 



Let A denote the lattice generated by the integrals of the basis of holo- 
morphic differentials around this basis of cycles in C. Then the manifold 
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C 6 /A is the Jacobian variety of C, denoted by J. Next, for k = 1,2,. 
define 21, the Abel map from the kth symmetric product Sym fc (C) to J. 



21 : Sym fe (C) -> J 



du-\ h 




du (mod A), (25) 



where the Pj are again points upon C. Denote the image of the fcth Abel 
map by W^ k \ and let [— l](ui, . . . , -u 6 ) = (— Ui, . . . , —u 6 ). We then define the 
kth standard theta subset (also referred to as the fcth stratum) by 



When k = 1 the Abel map gives a one dimensional image of the curve C. 
Since our mapping was given by a single integral with respect to one param- 
eter, it will make sense to rewrite this as an integral on the one- dimensional 
stratum, ©M. In addition to [4], [5] and [6], similar problems of inverting 
meromorphic differentials on lower dimensional strata of the Jacobian have 
been studied, in the case of hyperelliptic surfaces, in [14], [15], [16] and [17] 
for example. 

Let us also define a basis of second kind meromorphic differentials, dr, for 
the surface; these have their only pole at oo. These are determined modulo 
the space spanned by the du and can be expressed as 



A specific set was derived in [13] in order to construct Klein's explicit reali- 
sation of the fundamental differential of the second kind. This set was given 

as 

hi = -s 2 (l It 3 + 8t 2 /i4 + 5t/i 3 + 2/i 2 ) , h 2 = -s 2 (7t 2 + 4t/i 4 + /i 3 ) , 
h 3 = -2ts(3t 2 + 2tpL 4 + h 4 = -3ts 2 , h 5 = -2t 2 s, h 6 = -t 3 . 

We then define the period matrices u>', cu", rj' and rj" by 



Q[k] = W [k] y 




where 
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We can combine these into 

M 



u' uj" 
rj* rj" 



which satisfies the generalised Legendre relation, 

-!e A i\ jT in ( -1 6 



5 The Kleinian cr-function associated with C 

We will now define the multivariate cr-function associated with C, which is 
constructed from the ^-function, (see for example, [18]). This function is a 
generalisation of the classical Weierstrass elliptic cr-function, and as in the 
elliptic case, it can be used to construct Abelian functions on J. 

Definition 5.1. The Kleinian a-function associated with C is 

a{u) = a(u; M) = cexp ( - \ui{u')- l u T ) x 9[5] ((cu')- 1 ^ \ {u^lo") 
= cexp ( — iu?/^')" 1 -^) x ^ exp 

2vn j§(m + 5') T (cc/)-V'(m + 5') + (m + 5') 7 \{to')- l u T + 5") J 



Here c is a constant dependent upon the curve parameters, (/io, fii, fi2, f^3, A^)- 
The results of this paper are independent of this constant, so we do not discuss 

8' 
8" 

which gives the Riemann constant for C with respect to the base point oo and 
the period matrix [oj',uj"], (see [19] p23-24)- 



its value here. The matrix 5 = 



is the theta function characteristic 



We will evaluate the integrand using derivatives of cr(u), with respect to 
the variables u. We denote these cr-derivatives by adding subscripts. For 
example we write J-j^- as cr^. 

Lemma 5.2. We summarise the fundamental properties of the a-function in 
this lemma. Further details and proofs are available in [13]. For a detailed 
study of the general multivariate a-function, we refer the reader to [19]. 
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Given u G C 6 ; denote by u' and u" the unique elements in M 6 such 
that u = u'uj' + u"uj" . Let £ represent a point on the period lattice 

e = e u > + e'u" e a. 

For u, v G C 6 and lG A, define L(u, v) and x{£) as follows: 

L(u,v) = u T ( V 'v' + rf'v"), 

X (£) = exp [m(2(£' T 5" - £" T 5') + £' T £")] . 

Then, for all u G C 6 , £ G A £/ie function c(u) is quasi-periodic. 



a(u + £) = X {£) exp 



L(u+ £ -,£^ -a(u). (26) 



For 7 G «Sp(12, Z) we have 

a(u;-fM) = a(u;M). (27) 

o~(u) /ias zeroes of order 1 when u G and is non-zero elsewhere. 

In the case when all the curve parameters are set to zero, the function 
o~{u) is equal to a constant K, times the Schur-Weierstrass polynomial, 



- ^u\u% - ^sinful + Yfulul + | w 4 wi«3 - 2ulu 6 u 3 u 5 - u\u % 

- l u l u 3 u l - U A ul + ^u\ul - 3^4M6 M 4 ~ ^6 U l U 5 + lio8 U 6 U2 
+ \u\u2 + luluj - \u 4 ul + 39^68?4 2m 4 + U A U§u\u 2 + \u\ 

+ 2 U 5 U 3 U 2 + | U 5 2 U e 4 U2 + J2 U 6 5 U 2 U A - \ U 2 U§U 2 + \ U^U^U^ 2 

- | u 2 u§u^ — ^ m 5 4 m 4 -u 6 4 + U A U G Ui — U 5 2 Ui. (28) 

• The sigma function must have definite parity and weight. From SW45 
above we can conclude that o-{u) is odd with weight + 15. 

One of the key results in [13] was a Taylor series expansion for a(u) about 
the origin, used to derive relations between the Abelian functions associated 
with C . We will use this expansion in the evaluation of the integrand given 
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in Section 7. In [13] it was shown that this expansion could be constructed 
in the form 

O-(ll) = <J{U U « 2 , « 3 , «4, «5, «6) = Cl5(«) + Cig(«) H h Ci 5+4 „(w) + . . . 

where each C k was a finite, odd isobaric polynomial composed of sums of 
monomials in Ui of total weight +k, each multiplied by a monomial in fij of 
total weight 15 — k. From Lemma 5.2 we can conclude that C15 = SW^, 
while the other Cfc were found in turn, up to C59 by considering the possible 
terms, and ensuring the expansion satisfied known properties of a(u). (See 
[13] for full details of the construction, and [20] for a link to the expansion.) 

6 Deriving relations between the cr-function 
and its derivatives 

We will evaluate the integrand (21) as a function of cr-derivatives restricted 
to 0M. In order to achieve this we will need to derive equations that hold be- 
tween the various a-derivatives. In [13], sets of relations between the Abelian 
functions associated with C were calculated. However, these were the rela- 
tions that held everywhere on J, and do not give us sufficient information 
for the behaviour of o~(u) on the strata. To derive such relations we start by 
considering a theorem of Jorgenson [21]: 

Theorem 6.1. Let u E 0^ for some k < g. Then for a set of k points 
Pi = (ti, Si) on C we have 

u = ^2 du, 

i=i Jo ° 

and the following statement holds for vectors a, b of arbitrary constants. 

Z 9 j= i qgjfo) _ det [a\du{P^\ ■ ■ ■ \du(P k )\du(P k )^- k -^\ ■ ■ ■ \du(P k )^] 
E?=i^-(«) ~ det [b\du(P 1 )\---\du(P k )\du(P k )(9-^)\...\du(P k )W] 

Here, du^ denotes the column of ith derivatives of the holomorphic differ- 
entials du, and should be ignored if i < 1. 
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Below we state that the strata of C can be defined by the zeros of the o- 
function, and of its derivatives. The definition of 0® is a classical result (in 
Lemma 5.2) while the others can be derived from the theorem above. (See 
Appendix A for full details.) 



Q [5] 


= {u 


a(u) 


= 0} 








Q W 


= {u 


a(u) 


= cr 6 (u) 


= 0} 






Q [3] 


= {u 


a(u) 


= a 6 (u) 


= cr 5 (u) 


= 0} 




©[2] 


= {u 


a(u) 


= cr G (u) 


= ^(u) 


= (?a{u) 


= 0} 


QK 


= {u 


a(u) 


= <?%{u) 


= 05(11) 


= 0-4 (u) 


= 0-3 (u) 



(29) 



We use these defining relations to generate further relations between the 
cr-derivatives, holding on each stratum. We use a systematic method, imple- 
mented in Maple, to achieve this. 

Start with the relation cr(u) = valid for u G O' 5 '. Consider u as it 
descends to O' 4 '. We write u as u = u + it£ where u is an arbitrary point 
on and is a vector containing the series expansions, (24). We can 
calculate the Taylor series expansion in £ for a(u + u^) = as 

= a(u + us) = a(u) - a 6 (u)^ + \ [a m (u) - a 5 (u)]^ 2 + [\a m (u) 

- \(Ji{u) - |o- 666 (-U,)]£ 3 + [!<7 55 (w) + |(7 46 (w) - \cr 5(i6 (u) 

+ ^0"6666(^)]£ 4 + [j2 a 5666{u) + \cTAb(u) ~ ^(766666 (il) - \u Am (u) 

- i<j 5 5 6 (m) + ^ e (u)fi 4 ]e + o(e) 

Setting the coefficients of £ to zero gives us a set of relations for u e 6^: 

a 6 (u) = 
cr 66 (u) = a 5 (u) 

o-666(w) = 3<7 56 (u) - 2a 4 (u) (30) 
06666 (w) = 6o- 566 (w) - 8o- 46 (n) - 30-55(1*) 

0-66666 (u) = 10<75 666 (li) - 20cr 4m (u) - 15cr 556 (w) + 20<745(h) + 6^4 0"6 (u) 



If we calculate the expansion to a higher order of £ then more relations can be 
obtained. Note however, that the expansion for must first be calculated to 
a sufficiently high order first. We have calculated an expansion for a(u + u^) 
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up to 0(£ 29 ), using the weight properties of cr(u) to simplify the calculation. 
This expansion can be found on-line at [20]. 

The next step in this process will be to find the relations valid for u E G®. 
Since G^ C we can conclude that the relations (30) are valid here also. 
However, we can derive a larger set of relations for it E Q' 3 ^ by repeating the 
descent procedure for those relations that are valid on 6' 4 l 

We do not need to consider the relation u{u) = since that will only give 
us the same relations as above. Instead choose the second defining relation 
<T6 (it) = 0. We again write u = u + u$ where it£ is the vector of expansions 
as before and u is now an arbitrary point on 6®. We do not need to calculate 
the Taylor series expansion in £ for a 6 (u + it$) as before. Instead we can 
take the previous expansion and simply add 6 to each index: 

= a 6 (u) - (J 66 (w)£ + \ [a em (u) - a 56 (tt)]£ 2 + [\a m& (u) - \<Ju>{u) 
- !o"6666(w)]£ 3 + . . . 

Setting the coefficients of £ to zero gives us more relations valid for u E @' 3 l, 
starting with <7 66 (it) = 0. We can obtain further relations for u E 0^ by 
descending all of (30). We automate this process in Maple as follows: 

1. Take a relation valid for u E Q' 4 ' and expand as a Taylor series in £. 
To do this we replace each cr-derivative by the Taylor series expansion for 
a(u + u^), adding the relevent index to each a-derivative in the expansion. 

2. Set each coefficient with respect to £ to zero, and save the resulting equa- 
tions. 

3. Repeat steps 1 and 2 for all known relations valid for u E Q' 4 '. 

4. Use the set of equations we have obtained, to express the higher index 
ex-derivatives using lower-index derivatives. If we have a-derivatives with 
the same number of indices, solve for those with the higher indices first. 

Once we have finished this process we will have a set of relations valid for 
u E 0' 3 '. We can repeat the process by descending each of these to 
creating another set of relations which we can finally descend to 0W. We 
end up with a set of relations valid for u G 0' 1 ', some of which are contained 
in Appendix B, with the full set we have derived available online at [20]. 

The surprising result of these calculations was that on G)W, we have 
<Ji(u) = <7 2 (m) = along with the other first derivatives of a (it), concluded 
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to be zero using Theorem 6.1. These calculations were computationally much 
more difficult that in [6]. The latter stages were performed in parallel on a 
small cluster of machines, using Distributed Maple, (see [22] and [23]). 



7 Evaluating the integrand 

Recall our integrand, (21) 



2 . . 1 1 \ \ f (II 



<p(t)dt = K [Aj 1 + A 3 t + A 2 + -± + t2 { ^ 

Now, \{t) was given by a single integral with respect to one parameter, the 
point (t, s) on C. So we rewrite this as an integral on the one-dimensional 
stratum ©W of J, which we will parametrise by u±. We will then evaluate it 
using a-derivatives restricted to ©M. In [6] Jorgenson's Theorem was used 
to express t in terms of a-derivatives. However, if we solved (48) naively for 

u e 

eii] 

we would find 



t = 



which makes no sense given that <Ti(u) = for u E @M. Instead let us 
take equation (47) which was also derived from Jorgenson's Theorem (in 
Appendix A) and which holds for u on 0^. We consider what happens to 
this as u descends to ©W. We replace (£2,52) with the expansions in the 
parameter £ and replace the a-derivatives by their Taylor series in £. If we 
then take series expansion of this in £ and set £ = we find that for u G ©^ 
we have 

aia 23 (w) + a 2 cr u (u) _ ait - a 2 
M23(w) + b 2 a 34 (u) bit-b 2 

Solving this for t gives: 

< = "^M (3D 
a 34 (w) 

for u E @M. Therefore, using the basis of differentials (22) and equation (31) 
we can rewrite our integrand as: 

<p(t)dt = Kfatydt + <p 2 (t)dt] (32) 
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where 



<fi(t)dt — A 2 du 1 + A 3 du 2 + A^dui, (33) 

Mt)dt = ( (^r\) 2 - a 1 ^\) dui = Mu)d Ul . (34) 

y\a 23 (u)J cr 23 {u) J 

Thus <fi is a sum of holomorphic differentials on C , and ip 2 is a second kind 
meromorphic differential. As in the previous cases we will need to find a 
suitable function such that 

-^-*(u) = <p 2 (u), ^G0 [1] . (35) 

We will identify such a function ^f(u) as follows. First we must derive 
the expansions for (p 2 (u) at its poles. We will then find a function 
which has simple poles at the same points as the double poles of <p 2 (u), and 
which varies by at worst an additive constant as u moves round the a and 
/9-cycles of ©i. The function will be chosen so that -£^^(u), has the same 
expansion at the poles as (p 2 and is regular elsewhere. It then follows that 
the difference -£^^(u) — (p 2 (u) is holomorphic and Abelian; by Liouville's 
theorem we conclude that this difference is a constant. This constant may 
be evaluated at any convenient point. 



7.1 The expansion of ^{u) at the poles 

Recall that a(u) was an entire function, and so (f2(u) will have poles only 
when a 23 {u) = 0. Since we are working with u e 0' 1 ', by (31) this will 
occur at the points, one on each sheet, where t — 0. The cyclic symmetry [i], 
relating the different sheets of the curve, acts on (t,s) by [i] (t, s) i— > (t, is); 
hence, it will act on u as follows: 

u x im u 2 i m 2 u 3 h-> -u 3 , , 

U4 t— > IU4 u 5 I— > — m 5 -u 6 1— > — fu 6 . 

Let Wo be the Abel image of the point on the principal sheet where a 23 = 0. 
This is the point where t — and s = (fio) 1 ^. Then the full set of zeros of 
(723 are given by Wo,jv = [i^^o, N — 0,1, 2, 3. We will require the poles to 
match at all four of these points. 
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We need to find an expansion for ip 2 (u) at these points. To start, we 
consider u e 0^ and calculate the Taylor series of cr(u) around the point 
u = u Q = (wo,i,wo,2,«o,3,«o,4,wo,5,«o,6)- Writing Wi = (ui - u 0yi ) we have 

a(u) = a(u ) + [a 1 (u )w 1 + a 2 (u )w 2 + (T 3 w 3 (u ) + 0- 4 (u o )w 4 
+ a 5 (u )w 5 + a 6 (u )w 6 ] + [l<7 U (u )wf + a 12 (u )w 1 w 2 + o- n (u )wiw 3 
+ au(u )w 1 w 4: + 0-15(^0)^1^5 + Oi 6 (n )wiw 6 + |o- 22 (u )w2 + • • • 

We will have similar expansions around the other u ,n, and we can also use 
this expansion to easily compute the expansions for the o-derivatives, (by 
simply adding the relevant indices). Note, that since u 0jN are the points 
where t — 0, we can write their components as 



u o!i = / du h i = 1, . . . , 6. 

J 00 

evaluated on the sheet where s = [i] N (no) 1 ^ 4 . Therefore 

rt r 



ft pv ft 

[ui - v^J) = dui- du { = du u 

J co J 00 Jo 



evaluated on this sheet. Using (22), our basis of holomorphic differentials, 
we can find expansions for w^n, ■ ■ ■ w^n in the parameter t. 

1 i N _ 3 i 3N fi 1 +2 _ i 3N 8fi 2 fi -7fi 2 U 3 

Wl ' JV ~4/xo 3 / 4 32 /i 7 / 4 128 /i n /4 

_ %N (96/i 3 /i 2 - 168^i/i2/io + 77/ii 3 ) ^ 4 , ~ 

2048 J^m 1 +U ^> [6{) 

1 i N 2 1 i% - 3z^ 8/x 2 /i - 7/xi 2 4 ,~ 
2-^ ~~ sTTT374 "77^777771^ ~ 777^7 TTTfTi * + ° I* J 



/i 3 / 4 16 /i 7 / 4 512 /V 1/4 

1 j^, _ ifVi ,2 _ ^ 4^0-3 ^1^3 
^ 4 /i V2 t i 6 ^3/2 96 ^5/2 



Z 2jV 8/i 3 /i 2 - 12^1^2/iQ + 5/il 3 ^ / 5 x 

okr ,,7/2 1 + u V ) 



256 ^q 7 / 2 

3 



(38) 



21 



^"8^ "24^ " 128 ^5/2 ' +°^) 

1 ^\ _ _ 8/i2/iQ ~ 5/i! 2 ^ 

W6 ' JV "4 // V4 t 32 384 /.o 9 / 4 

_ » 3JV (32/i 3 /iQ 2 -40/x 1 /x 2 /x + 15/i 1 3 ) 44 , 5 . 
2048 ^ + 1 > 

Note that all these expansions are given for the general sheet since we need 
to check the behaviour at all the poles. We can move between the sheets by 
selecting the appropriate value of N. We can invert (37) on the iV-th sheet to 
give an expansion for t in tu^jv, allowing us to use w^n as a local parameter 
near u ,at- 

t = 4 z 3A Vo 3/ V,;v + 6^ N ^ /2 wl N + O (wl N ) . 

We start by substituting for t to give the expansions of W2,n, ■ ■ ■ , we,N with 
respect to w^n- 

w 2 , N = 2i 3N fi 3 ^wl N + O (wl N ) 

; (39) 

■2N 1/2 1/4 -5N 2 , ^ ( 3 \ 

W 6 , N = I /J, ' w hN + n ' Hxl W 1N + O {w 1N ) . 

We use these in turn to give the ex- derivative expansions at uo,jv as series in 
Wi ; n. For example we have, 

cr 23 (w) = cr 23 (wo,jv) + (« 2A VoO"236 + ^123 + i^o ^) ( w o,n)«Ji,jv 
+ O Kn) ■ 

We substitute these into (34) to obtain an expansion of ^(u) at u — u 0jN , 
as a series in w 1 jv- 



0"34 



,2 



1/4 -2^, 1/2, 



W 1,AT V^ V //q i(T 233 + ^ J >0 "0-236 + ^123 
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where C(u 0jN ) is a polynomial in the a-derivatives, which we need to evaluate 
to ensure that </? 2 (ii has zero residue. 

In the previous section we derived a set of relations for u E O^, but 
these are not sufficient to simplify C(wo,jv)- We need to generate a further 
set of relations which are valid only at u — uo,jv- We do this using a 
similar approach to the previous section. We take a relation valid on 0W 
and calculate its expansion around u = wo,at as a series in Wi^n, using the 
series derived above. We then set to zero the coefficients of u>i 5 at. We do 
this for each relation valid on @M and obtain a set of equations between 
cr-derivatives at the points u = wo,at- The first few such relations are given 
in Appendix C, with a fuller set available online at [20]. 

If we substitute these into the expansion of <^2( u ) we obtain 



ip 2 {u) = 



■2N 



16 fl* 



W 



+ 



1,N 



;N 



— (4/^1 - 3/^) 



Win 



+ 0(w%) (40) 



Recall equation (18) which stated that ip(p) has zero residue at p = oo on all 
sheets. Since residues are invariant under conformal maps, we can conclude 
that ¥2(11) must also have zero residue, and so the constant A\ must be equal 
to 



A, --^ 
4 no 



giving ip 2 (u) = 



•2 AT 



16 »r 



i 



w 



+ 0(w° hN ). (41) 



1,JV 



7.2 Finding a suitable function ty(u) 

We need to derive a function \J/(it) such that the Laurent expansion of 
has the same principal part at the poles as <p2(u), so we will restrict our 
search to linear expressions in a-derivatives, divided by a 2 3(u). For these 
functions we will derive expansions in w^n at u = u ,jv using the techniques 
described above. Let us take the function 



*(«) 



Ecrdu 
n 

1<!<6 



<Tij[U) 



a ^ U ) l<^<6 a ^ 



E 



Vijk 



l<i<j<k<6 



where the rjij and rjijk are undetermined constants. (Through trial and error 
we found that 3-index cr-derivatives were required in the numerator). Now, 
since we are working on we will find that many of these a-derivatives 
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are equal to zero, or can be expressed as a linear combination of other such 
functions using the equations in Appendix B. Let us set the coefficients of 
these functions to zero, leaving us with 

Vj/( W ) = [rj 22 (T 22 + + T/lllOlll + ??122 0122 + ??1230"123 + ^1340"l34 

+ ??2220"222 + 7/223C223 + ??2240224 + ^225°"225 + ^226°"226 + ??233 0233 
+ ^234°234 + V235° "235 + ??2360236 + %34°"334 + ^344°"344 + ^345 0345 

+ 773460-346] («) 



0-23(1*)' 

We emphasise that we need to work with the total, not the partial, derivative 
of ^/(u) with respect to Mi; in practice the other Ui are expressed in terms of 
w± t N in the vicinity of n. ,iv so there is no ambiguity. Note from (22) that 



-2- =1 
Therefore 



' dui ' 



d 
duz 



dui 



d 
8u4 



t 2 — 

dui 



d 

(9«5 



St 



dui 



d 
due 



2 d 
dui ' 



D d 

1 dui 

d 



d 

_ g- 23 (tt) g 
(9-Ui a u (u)du 2 



t 



d 



+ s- 



d 



d 



du 2 du 
d 



+ st 







d 



+ s 







du 5 du e 



— s- 



023 9 



023 (Wj 



034 («)/ <9m 4 0-34 <9« 5 



We can now evaluate ^-^(u) as a sum of quotients of cr-derivatives. For 
example 



+ 



0236 
023 
°"23 er 2346 
3 2 4 



01236 02360123 



023 (T 

02360234 
3 2 4 



2 

23 



02236 . 0236 0223 . 02336 0236 0233 
1 h S S 7, 



034 



023034 



023 



'23 



02356 . 0236 0235 . 2 02366 

S k h s h s 



'34 



023034 



023 



2 

,2 0236 



'23 



Now let us consider the expansion of Di(^(u)) at u — u 0jN . We generate 
series expansions in tu^jv for the relevant a-derivatives using the method 
described in the previous subsection. We can use the relations in Appendix 
B and C to simplify these expansions, and so obtain a series in wi t N for 
Di (*(«)). We find, 



U=U N 



N(u , N ) 

022 (WO, Ar) 



it; 



1,7V 



+ 0(wl 



N)i 
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where N(u , N ) is a linear polynomial in {a 22 , cr 122 , a 222 , <t 223 , <t 224 , <t 225 , <j 226 }. 
This set of a-derivatives can be used to express all other 2 and 3-index o- 
derivatives when u = Uq^n, (as in Appendix C). We find the coefficients 
of N(u ,n) with respect to each of these seven a- derivatives and determine 
conditions on the constants that set all the coefficients, except that 

of cr 22 , to zero. We then obtain further conditions by ensuring the expansion 
we are left with (now independent of any cr-derivatives) is equal to (40) on 
the four sheets. Imposing these conditions on \I/(it) leaves us with 



C22 , fin <7 2 35 8?7 22 /io — 1 0"236 , C"334 

V22 h 77111 h 2T/334 1 V ??334 

023 C"23 C"23 4/iQ 0"23 ^"23 



U 



Note from Appendix C that the terms containing cr m , cr 235 , 0334 all vanish 
at the points u = uo,n and so have no effect on the expansion here. Let us 
discard these to leave, 

lTr / \ cr 22 (u) 1 (8r/ 22/ u - 1) cr 236 (ti) 
= 7/22 7-t + 7- 



'cr 23 (u) 4 /i (723 (w) 

We now have two functions, v 9 2(^) and Di(\E'(n)), which both have poles 
at uq^n- We have derived expansions at these points, given in the local 
parameter u^jv, and ensured that they match. However, we should also 
explicitly check what happens at the point u = 0, since tu^jv is not a suitable 
local parameter here. We can instead use the Taylor series expansion of cr(u) 
presented in [13] and described in Section 5. We differentiate this to give 
expansions for the cr-derivatives and then, since we are at the origin, replace 
the variables u±, . . . , u 6 with their expansions, (24), in the local parameter £. 

Now, the sigma expansion was given as a sum of polynomials with in- 
creasing weight in u and hence the expansions will have increasing order of 
£. Since the functions we consider all contain ratios of cr-derivatives we will 
only need the leading terms from each expansion, in order to check regularity. 
Hence we only require a minimum amount of the sigma expansion, sufficient 
to give non-zero expansions for the derivatives we consider. We find that for 
the functions used here, we can truncate the expansion after C35. 

Substituting these expansions into (/^(tt), we find 



lim v? 2 (n) = lim 



3/il 4+ 8 + 36 

4/V 



0. (42) 



So v ? 2(ii) is regular at the origin, and hence we must ensure that fy(u) is as 
well. Upon substitution into ^f(u), we find that we must set r/ 22 = for the 
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expansion to be regular. This leaves us with 




1 1 <T23(i(u) 

4/i <t 23 (u) 



with 




1 (-8/x 2 + 3/i 4 /i 3 ) 
176 fi 



£ n + 0(£ 15 ) =0. (43) 



Now all that remains is to check the periodicity properties of the functions. 
Recall equation (26) which gave the quasi-periodicity property of a(u). We 
can differentiate and use the relations in Appendix B to show that 



Hence, from equation (34) we can see that <p 2 (u + £) — (p 2 (u) and so is 
Abelian. Similarly we can see that 



Hence Di(^(u)) is Abelian, and we may now write our integrand ip(t)dt as 
Lp 2 (u)dui + A 2 du 1 + A 3 du 2 + A A du A = + B T du, (44) 

for some vector of constants B T = (B 1: B 2 , B 3 , B 4 , B 5 , B 6 ). 

7.3 Evaluating the vector B 

We can evaluate the vector B by considering the integral of equation (44) at 
the point u = 0. We can use the expansions in (23), (42) and (43) to obtain 



a 34 (u + £) 
a 236 (u + £) 



X (£)exp(L(u + l,£))a 23 (u), 
X (£)exp(L(u + l,£))a u (u), 

X (£) exp(L(u + f , £)) ■ [JLL(u + a 23 + a 236 ] . 



^(u + £) 



1 1 a 236 (u + £) 
4/i a 23 (u + £) 



*(u) 



1 1 d 

4 no du 6 
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the following series. 

+ 7^ (3/X 4 /X ^4 — A*3 — 3/i4/Uo5 4 - 4A 3 /i + 



(32A 2/ u - 325^0 + 16^2 - 6^4/13 - 24fi A 4 fi 3 + 2l^A A n\ 



352 /i 

+ 24/i 5 4 /i3 - 21/i 5 4 /i4 ! - 24/i 4 /ioA 3 + 24/i 4 /i 5 2 )£ 11 + 0(£ 13 ) 
Setting each coefficient of £ to zero, we find 

1 (/i 2 + 2A 2 ^) 1 (4A 3 /i + /i 3 ) R „ 

2 no 4 /i 

£? 4 = A 4 , B 5 = 0, £? 6 = 0. 

8 Obtaining an explicit formula for A(p) 

We now use the results of Section 7 to derive an explicit formula for the 
mapping X(p). Start by applying the change of coordinates given in (19) to 
X(p) as given in (15). 



Hp) = v + £w) - iW = (ft (v® - 1) 







+ ^ 8 ~ P (k[A4 a + A 3 t 3 + A 2 t 2 + Ait + , 

where the constants A±, A 2 , A 3 , A4 and K were defined by equation(21). Note 
from (31) that 

P = Ps+ a -^r\, UG0N. (45) 

023(«) 
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So let us take u e 0' 1 ', and use (45) and the evaluation of (p(t)dt from the 
previous section to write A(p) as 



A(p) = [P8 + 



cr 34 (w) 



+ 4 



023 («) 

1 (4A 3 w + /is) 



Wo 

Integrating, gives 

0-34(1*) 



du 2 + A^dui 



1 

T 2 

+ K 



dt + K 



f P8~P 

Jo 



d 



1 (// 2 + 2A 2< U ) 

2 u 

1 1 o 236 (u) 



dwi V 4/i (72 3 (w) 



A(p) =[P8 + 



+ 



0"23(w) 

1 (4A 3 u + "3 



"o" 34 (w)~ 




.023 





4 /i 

= P 8 + X 



W 2 + 



1 (/x 2 + 2A 2 fx ) 

2 m 

1 1 O 230 (li) 



-Ml 



4u o 23 (ii) 



1 o- 236 (u) u 2 + 2A 2 u 4A 3 /i + u 3 

H Mi H u 2 + A 4 u A 



4/i Cr 23 («) 



2/i 



4«o 



+ C. 



for some constant C. We can determine C by ensuring that the following 
condition on the mapping is satisfied. 

lim A(p) =p + ( - 

Note from (45) that p — > 00 implies <7 23 (it) — ■> and therefore w — > mo,jv- 

o 34 (w) 



lim [A(p) — p] = lim 



c- 



+ K - 



23 (Wj 

1 1 o 236 (u) « 2 + 2A 2 /i , 4A 3 « + «s 



4 w o- 23 (u) 



+ 



2/i 



-Mi + 



4w 



U 2 + A4U4 



Let us ensure that the condition is met on the sheet of the surface C asso- 
ciated with limf^ ( s ) — We can write this as a series expansion in the 
local parameter w±, (as described in Section 7.1). Recall that u,i = Wi + u ^, 
and use the expansions (39) and the existing expansions for the o-derivatives 
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to obtain 



lim [A(p) -p] 

p—fOO 



11 K 1 

ij ~ 16" 



A*6 



1 



8 no 



1 /ii (^2 + 2^2/10) (4A 3/ u + /i 3 ) 
+ 32^ + 2^ ^ + 4^ M0 ' 2 + 

Therefore, we must set the constants of integration, C, to be 



1 1 a 22(i (u ) 
4// cr 22 (uo) 



3 /i 



C = -— + K\ +- 



1 1 (7226 (Uq) 1 A*l 1 (// 2 + 2A 2 /X ) 



8 no 



4// (T 22 (uo) 32^/4 2 



-«o,i 



1 (4A 3/ u + At 3 ) 



U ,2 - A 4 U ,4 



4 /i 

giving us the following explicit formula for the mapping X(p). 



+ 



3 h 
8^ 

fi 2 + 2A 2 /j, 



/io 



1 1 



^236 (^) _ q-226(^o) 
CT 23 (w) Cr 22 (u ) 



1 /il 



2/i 



4A 3 /x + /i 3 

(Ml - Mi i0 ) H [U 2 - U 2 ,0) + A 4 {u A - U 4 fl) 

4/i 



where u e 6^ and u is the point on the principle sheet of the surface C 
where t — 0. 
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A Deriving defining relations for the strata 
from Jorgenson's Theorem 

Consider Theorem 6.1 in the case when k = 5. Then for arbitrary a, b 



E 7 -=i^j(«) 



det [a 


du(Pi) 




du(P 5 )] 


det [b du(Pi) 




du(P 5 )] 



Now, as u e 0® approaches 0^ we have the point P5 = (£5, S5) approaching 
00. We can use the local coordinate £ here and hence replace the final column 
of the determinants by the expansions (23). When u arrives at 0' 4 ' we will 
have £ = and hence the determinant in the numerator becomes 



a 2 
a 3 

0,4 

a 6 



dt ± 
ti dti 



t'fdh 
4 S 3 
tisidti 

s'fdti 
~4&f 



din 
S4ai4 



4s? 



Is 



^4' 



46?f 



4sf 
















a\ 


1 • 


1 









a 2 


h ■ 


• t 4 





dti 


dt^ 


a 3 


si ■ 


s 4 





3sf ' 


'3sJ 


0,4 


ti ■ 


t 2 









a 5 


t 1 s 1 ■ 


■ £484 









a 6 


si ■ 


s 4 


-1 



The determinant in the denominator will be identical, except with the entries 
of a replaced by the entries of b. Hence the factored terms will cancel, 
leaving us with the simpler determinants. It is clear from the final column, 
that when we expand the determinants the resulting quotient of polynomials 
will not vary with the arbitrary constant clq. Hence we must conclude that 
for u E 0W, <7&(u) = 0. The same conclusion could have been drawn from 
considering 6 6 . 

eW = {u I a{u) = <7 6 (u) = 0}. 
We repeat this process by considering Theorem 6.1 in the case when k = 4. 



1 a j (T A u ) 



det [a 


du(Pi) 




du(P 4 ) 


du(P 4 )( 1 )] 


det [b du(Pi) | ••• 


du(P 4 ) 


dn(P 4 )( 1 )] 



This time we consider u descending to 0' 3 ', by letting the fourth point move 
towards infinity. The penultimate column in each determinant can be given 
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with the expansions (23) as before. For the final column we will need to 
determine the derivative of these expansions: 



d 2 Ui 

W 
d 2 u 2 

d 2 u 3 



-10£ 9 + O(£ 10 ) 
-6£ 5 + 0(£ 6 ) 
-5£ 4 + 0(£ 5 ) 



^ = -i + lM 4 +o(a 



(46) 



cPuq 



M 3 + o(a- 



When u arrives at 6' 3 1 we will have £ = 0. Our determinants will again 
factor and cancel to leave the numerator as 





1 • 


■ 1 








a 2 


h ■ 


• t 3 








03 


si • 


■ s 3 








04 


t\ ■ 


t 2 








a 5 




■ hS3 





-1 


a 6 


si ■ 


S 2 


-1 






with the denominator identical with b instead of a. From the final two 
columns it is clear that the resulting quotient of polynomials will not vary 
with the arbitrary constant a 6 and a 5 . Hence we must conclude that 

= {u I a(u) = a 6 (u) = a 5 {u) = 0}. 

We repeat the procedure once more for k — 3. We let u descend to 0^1 and 
use the previous expansions along with the derivatives of (46) for the final 
three columns. We let £ = 0, cancel the common factors and expand the 
determinants to reduce the statement to 



(47) 



Yfj=i a j< J j( u ) _ aiU s 2 - oi£j t 2 + a 2 5j - Q2gg - azh + a 3 t 2 
Ej=i bj(Tj(u) hh s 2 - hs! t z + b 2 s 1 - b 2 s 2 - ht t + b 3 t 2 

for u E 0l 2 l. We therefore conclude that 

0[2] = [ u I a(u) = a 6 (u) = a 5 (u) = a 4 (u) = 0}. 

Finally we consider Theorem 6.1 in the case when k = 2. Here, when we 
let u descend to 0M, we find that the statement of the Theorem involves 



31 



singular matrices, (resulting from the final set of series for the derivatives all 
equalling zero when £ = 0), and hence gives us no information. 

Instead we can consider equation (47) which held for u e 0' 2 '. Let u 
descend to 9^ here, by using the expansions in £ for (£ 2 , s 2 ). We find that 



£?=i a j°j(«) Mi 



a 2 



h ( \ ~ bt b + °® (48) 

and so for u e we can see there is no dependence on 03,04,05 or a®. 
Hence 

©W = { w J er(it) = (T6( w ) = &b{. u ) — &4:(u) = cr 3 (w) = 0}. 
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B Relations between er-derivatives on 

[i] 



The following list of equations are valid for u £ @M. This set contains all 
those relations we have obtained that express n-index u-functions for n < 4. 
A larger set that includes relations for n > 4 is available online at [20]. 



^ 1 


: 


(Tii = 0, 


024 = 





<j 44 = 






(To = 


: 


^12 = 0, 


025 = 




(J A r- = 
^ 40 








= 


0"13 = 0, 


026 = 





(T^e = 






(J A = 


= 


014 — — 2 <T22 ' 


033 = 





CTt-c: = 






^ 


= 


015 = -023, 


035 = 





CTt-o = 
00 






(Jfi = 
o 


= 


0-16 = 0, 


036 = 





(Jo R = 
DO 






0112 




^0°34 + ^1023 




0156 


= —0236 


0445 


= 


0113 









0166 


= 023 


0446 


= 


0114 




-0122 + /^1034 + /^2023 


0244 


= 023 + /^4034 


0455 


= 


0115 




— 2(7123 




0245 


— ~2 0344 


0456 


= 


0116 







1/^3023 


0246 


= 


0466 


= 


0124 




— 10222 + |^2034 + 


0255 


= -2o- 345 


0555 


= 


0125 




— 2 0223 ~~ ^i 34 




0256 


= —0346 


0556 


= 


0126 









0266 


= 034 


0566 


= 


0133 









0333 


= 


0666 


= 


0135 




1 _ 

"2 0233 




0335 


= 






0136 









0336 


= -2(T 2 3 






0144 




-0224 + At4023 + At3034 


0355 


= 






0145 




— 0234 — 2 CT225 




0356 


= -034 






0146 




— 2 (T226 




0366 


= 






0155 




— 0334 — 2(T235 




0444 


= 3(734 







01136 = — /^O034 

<7n44 = /^30223 _ | ff 2222 + 2(/! 4 <7i23 + /^20234 + /^30134 ~ 01224) + /^10344 

01145 — /^20235 — 2(7i234 — 01225 _ ^02223 + 2^3 0233 + /^20334 + ^10345 

01146 — — 01226 + ^2 022 + ^2^236 + ^10346 

01155 — — 2(7i334 — 4(7i235 — 02233 + ^2^22 

01156 — —^1034 — 2(7i236 
01166 — 2<7i23 

(71244 = —§02224 + 0123 + fA*40223 + /^40134 + |/^20344 + /^30234 
01245 — 4 ^40233 ~~ q 02225 ~~ o" " 2234 ~~ 2°" 1344 2^ i3 °" 334 
+ 1/^2 0345 + ^30235 
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0"1246 = 2/ i 3°"22 — g°"2226 + fA* 3 * 7 ^ + I^ 20 " 346 

°"1255 — ^3^22 — 0"2235 ~~ 27 1345 — 02334 

C"l256 — — 2 " 2236 _ 2^ 2 °" 34 ~~ °"1346 

^1266 — 2 a223 + a 134 

Cl333 = -6/ioO"34 + 2^i<723 

4 4 1 

°"l335 — 3/^2^23 — 3/^1^34 — 3^2333 

C1336 — — 20123 

°"l355 = "02335 — o°"3334 _ ^2^34 + ^^3 a 23 

1 1 
°"l356 — — t" " 223 _ °" 134 _ 2 " 2336 

°"l366 — 2 0233 

01444 = §0223 + 30134 — §02244 + 3/i40234 + §/i30"344 

01445 = /i3<J345 + \<J233 — 02245 — 02344 + /^40334 + ^4^235 



01446 




114&22 — 02246 + A*30"346 + /M0236 






01455 




— 1^3344 — 2(T2345 — ^02255 + ^4022 






01456 




— 2°"2256 — 2"/ 130 " 34 ~~ °" 234 6 








01466 




^"234 — 2 " 2266 








01555 




—302355 + 8/14023 — 803345 








01556 




—2(7234 — 0"3346 — 2o"2356 








01566 




0334 — 02366 + 0235 


03556 




— 27345 


01666 




022 + 3(7236 


03566 




— 27346 


02444 




3o"234 + §A t 40"344 


03666 




0"34 


02445 




0334 + °235 — 3 1 03444 + ^4^345 


04444 




6(7344 


02446 




022 + 0236 + A*4C346 


04445 




37345 


02455 




022 — 03445 


CT4446 




37 346 


02456 




— 2^4^34 — 2°"3446 


04455 







02466 




1 _ 

2 a 344 


04456 




-034 


02555 




10(T23 + 2/i4<T 34 — 3CT3455 


04466 







02556 




— 0344 - 2(73456 


04556 







02666 




3<7 3 4 6 


04566 







03333 







04666 







03335 







05555 







03336 




— 3(7233 


05556 







03355 







04555 




4734 


03356 




— 2<7 33 4 — 27 2 35 


05566 







03366 




-27 22 - 4(7236 


05666 







03466 




0345 — 02566 


06666 







03555 
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C Relations between a- derivatives at ito,iv 



The following list of equations were not valid in general for u £ 0W, but are 
true at u = u 0jN . This set contains all those relations we have obtained that 
express n-index cr-functions for n < 4. A larger set that includes relations 
for n > 4 is available online at [20]. 



1 022 



0m 
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0113 
0123 
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» 022/ii 
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■N 1/4 
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1 (X222 1 ^2022 

1 0223 



4* 2Ar ^ 
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"-022 
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1 022/^3 0224 
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1 0226 
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01112 = 
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2 ^ttf 

3 • 9 AT 

-6/i O 2 022« 
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A 1 ! 0"22 



3., 5/4-. -AT 
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— I V/io/i2022 



3 3/Xi(7i22 , 1 3 

" 2 (Tll22 + 2"^ + 2 /i04 ° r222Z 



3 AT 



7V 



2V 



2 022/^0 



1„'3AT 



^0 4 0223 - 



/iO« 2W 0122 
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3 « 2JV /x 2 cr 22 /ii 3 . 2N 
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1 A t lO"222 ^ 1 /i2°"l22 



01126 
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1 
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3Af 



6 
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3 0"2233 5 ^20"22 
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